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ABSTRACT 

We analyze the Suzaku XIS data of the central region of supernova remnant G332.5-5.6. 
The X-ray data are well described by a single non-equilibrium ionization thermal model, 
vnei, with an absorbing hydrogen column density of 1.4lg ^ x 10 21 cm -2 . The plasma is 
characterized by an electron temperature of 0.49lgQg keV with subsolar abundances for O 
(0.58lg 05 solar value) and Fe (0.72 +q'q| solar value) and slightly overabundance for Mg 
(1.23lg '^4 solar value). It seems that the central X-ray emission originates from projection 
effect or evaporation of residual clouds inside G332.5-5.6. We estimate a distance of 3.0 ± 
0.8 kpc for G332.5-5.6 based on the extinction-distance relation. G332.5-5.6 has an age of 7 
- 9 kyr. 
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1 INTRODUCTION 

Supernova remnant (SNR) G332.5 -5.6 was identified as an SNR 
candidate by iDuncan et alj (1 19971) based on t he morphology at 
4.85 G Hz and the lack of 6 0 //m emission. iRevnoso & Greenl 
J2007I) and lStupar et al.l J2007t) both analyzed the 1.4 and 2.4 GHz 
radio continuum data from Australia Telescope Compact Array 
and obtained an average radio s pectral index of a = -0.7 ± 0.2 
and a = -0.6 ±0.1 respectively. IStupar et al. 1 120071) analyzed the 
optical spectra of G332.5-5.6 and gave average ratios of [Nll]/Ha 
= 2.42, [Sll]/Ha = 2.10 and [Sll]6717/6731 = 1.23 which are all 
well within the typical values of SNRs. Due to these evidences, the 
authors confirmed G332.5-5.6 as an SNR. 

SNR G332.5-5.6 has an angular size of about 30' at radio 
wavebands with three main structures: two straight and nearly par¬ 
allel outer filaments at northeast and southwest, an extended radio 
emission region at center with small angular i nternal substructures 
dRevnoso & Greenl 120071 : IStupar et al.l i2007l) . G332.5-5.6 shows 
highly linear polarization of 0.3-0.5 at 2.4 GHz. IRevnoso & Greenl 
j2007l) mentioned that the high polarization might be artificial, 
because the Q/U intensity images might have more small scale 
structures than the total intensity image and the lack of short 
baselines excludes more total intensity emission than Q/U intensity 
emission. ROSAT PSPC observation reveals strong X-ray emission 
from the center of the remnant, which has similar morphology 
with the central radio emission. Besides a small bright point-like 
source in the northeastern filament, the extended X -ray emission in 
the tw o outer filaments is weak (see Figure 7 of lRevnoso & Greenl 
120071 and Figure 17 of IStupar et~aJ] |2007l) . The central X-ray 


emission can be well fitted by the thermal Raymond-Smith 
model with a temperature of about 0.25 keV l lStupar et al.ll2007l) . 
Neither a radio nor a X-ray pulsar is found in the vicinity of 
SNR G332.5-5.6. It seems that the central X-ray emission is not 
from a nonthermal pulsar wind n ebular (PWN) but from cloud 
evaporation or thermal conduction (iRevnoso & Greenll2007l) . 


Ha filaments are detected mainly in th e cen tral region and 
match well with the radio structure Jstupar et al.l2007t ). The spectra 
diagnosis of [Sll]6717/6731 suggests an average electron density 
of about 240 cm -3 . The H/3/Ha flux ratio indicates a reddening 
of E(B — V) = 0.27 towards G332.5-5.6. These facts indicate the 
blast at the center of SNR G332.5-5.6 is propagated into a cool 
and dense environment. Distance to the remnant is estimated in 
the range of 3 to 4 kpc based on Hi a bsorption, reddening and 
the surface brightness-diameter relation iRevnoso & Greenll2007L 
IStupar et ~al1|2007h . 


Our knowledge about G332.5-5.6 is still poor. In this paper 
we present the Suzaku X-ray Imaging Spectrometer (XIS) obser¬ 
vations of SNR G332.5-5.6 to determine its basic parameters, e.g. 
shock velocity, age, evolutionary state and elemental abundances. 
We also re-estimate its distance and discuss the origin of its central 
X-ray emission. 
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2 DATA AND DATA REDUCTION 

The central region of SNR G332.5-5.6 was observed by the XIS 
on board of Suzaku , starting on 2007 August 18 at 16:39:41 
and ending on 2007 August 19 at 20:57:24 (observation ID 
502066010, PI: Reynold S.). The pointing direction was (Ra, Dec) 
= (16h42m56.04s, -54d30m48.6s). Only the CCD camera XIS0, 
1 and 3 were used during the observation with clocking mode of 
8s exposure per frame and editing mode of 3 x 3. Each CCD 
has field of view of 17.8' x 17.8 / and angular resolution of about 2'. 

Data reduction is performed using HEAsoft version 
6.15. We reprocess and screen the data with the command 
Aepipeline by applying the up-to-date calibration file and 
then obtain an effective exposure time of 70.15 ks, 70.15 ks 
and 70.16 ks for XIS0, XIS1 and XIS3 respectively. XSELECT 
version 2.4c is employed to extract images and spectra. We used 
Xisnxbgen to generate the XIS non-X-ray background (NXB) 
image and spectra based on the night Earth data. The rmf and ar f 
files are created by Xisrmfgen and Xissimarfgen respec¬ 
tively. To fit the spectra, we use the software XSPEC version 12.8.1. 

3 RESULTS 
3.1 X-ray Image 

Figure[T|.shows the XIS0, 1 and 3 combined images of G332.5-5.6 
in the energy band of 0.5-1.0 keV (a), 1.0-2.0 keV (b), 2.0-3.0 keV 
(c) and 2.0-5.0 keV (d). Figure [2] displays the 843 MHz continuum 
image of G 332.5-5.6 (data from the Sydney University Molonglo 
Sky Survey. iBock et al]| 1999h with the X-ray bright region labeled 
by an ellipse. X-ray emission of G332.5-5.6 is distinct in the 
0.5-1.0 keV and 1.0-2.0 keV bands with their morphology very 
similar to the extended central radio emission. There is nearly no 
emission above 2.0 keV. The central X-ray emission has two bright 
substructures centered on (Ra, Dec) = (16h43ml5s, -54d33ml7s) 
and (16h42m35s, -54d29m33s) respectively. The northwestern part 
of G332.5-5.6 marked by the white curves in Figures Q]c) and (d) 
is contaminated by stray light from a low mass X-ray binary 4U 
1636-536 which is bright and about 50' away from the G332.5-5.6. 
Further spectra analysis has avoided that region. 


3.2 X-ray spectra 

Figure QJa) displays the region where we extract the spectra 
of the SNR (green dashed polygon, hereafter on-SNR). The 
on-SNR is outside the region contaminated by stray light. The 
background spectra are extracted from a region outside the SNR 
(hereafter off-SNR) and within the field of view of observation 
ID 502067010. The observation ID 502067010 points to (Ra, 
Dec) = (16h44m22s, -54d22m42s) and covers part of the north¬ 
west filament of G332.5-5.6. To increase the statistical quality, 
the front-side illuminated CCDs XIS0 and XIS3 spectra are 
merged. Then the SNR and the background spectra are binned to a 
minimum of 25 counts per bin to achieve a signal to noise ratio of 5. 

Since we do not care for the real ingredient of off-SNR, 
the NXB - subtracted XIS0+XIS3 and XIS1 off-SNR spectra are 
fitted phenomenologically and simultaneously with the model 
phabs* (vnei+powerlaw) within 0.4 keV to 4.0 keV. The 



Figure 1. The Suzaku XIS X-ray image of G332.5-5.6 in the energy band 
of 0.5-1.0 keV (a), 1.0-2.0 keV (b), 2.0-3.0 keV (c) and 2.0-5.0 keV (d). 
These images are smoothed with a Gaussian distribution of a = 24”. The 
green dashed polygons in (a) shows the on-SNR region. The white curves 
in (c) and (d) represent the regions affected by stray light. The ellipse in (b) 
occupies the same area as the ellipse in Figure [2] 



Figure 2. The SUMSS 843 MHz image of G332.5-5.6 with an ellipse indi¬ 
cates the central X-ray emission region. 

fitting is showed in Figure[3]with the reduced \ 2 (dof) = 1.15 (198). 
We amplify the off-SNR model by area ratio between on-SNR and 
off-SNR and treat it as a fixed component in on-SNR spectra fitting. 

Figure [4] shows the NXB-subtracted spectra of on-SNR. Line 
features (stronger than those in the off-SNR spectra) from O VII 
(0.57 keV), O VIII (0.65 keV), O VIIl/Fe XVII (0.81 keV), Ne 
IX (0.92 keV)and Mg XII (1.35 keV) are clearly showed in the 
spectra, indicating thermal emissions from the SNR. We first try 
to fit the spectra with single collisional ionization equilibrium 
(CIE) models, vapec and vmekal. After several tests, we find 
that none of them can give acceptable fitting to the spectra (Both 
CIE models have reduced \ 2 > 1-43 and can not describe the 
spectra well in energy lower than 0.6 keV.). Thus we then turn 









Suzaku observation and distance of G332.5-5.6 3 



Figure 3. The NXB-subtracted XIS0+XIS3 (black plus) and XIS1 (red 
plus) spectra of off-SNR. The solid black and red lines show the best-fit 
result. 


Figure 4. The NXB-subtracted XIS0+XIS3 (black plus) and XIS1 (red 
plus) spectra of on-SNR. The solid black and red lines show the best-fit 
result using vnei by changing the column density of hydrogen. 


to the non-equilibrium ionization (NEI) model, vnei, which is 
characterized by a constant temperature T and single ionization 
timescale r with flexible elemental abundances. The ionization 
timescale is defined as r = n e t, where n e is the post-shock electron 
density and t is the time since the passage of the shock. We set the 
vnei version to 3.0 which uses the beta release of AtomDB 3.0 to 
calculate the spectra. 

At the beginning of the vnei fitting, both the plasma 
temperature T and ionization timescale r are varied. The absorbing 
column density of hydrogen nuclei Nh is fixed to 1. 5 x 10 21 cm' 3 
based on the relation Nh = 1-79 x 10 21 Ay liPredehl & Schmitt 


fitting is poor with reduced x 2 (dof) = 1.63 (536). Then we try to 
change the abundances of some of the elements in the spectra. 
Changing the abundance of oxygen can improve the fitting evident 
with a new x 2 (dof) of 1.46 (535). The fitting can be further 
improved by changing the abundance of iron with x 2 (dof) = 1.30 
(534). Varying the abundances of magnesium improves the fitting 
with x 2 (dof) of 1.28 (533). Changing the abundance of other 
elements, such as neon, gives little or no improvement to this 
fitting, therefore their abundances remains fixed to solar value. We 
also try to free the column density of hydrogen, then the x 2 (dof) 
is 1.29 (532). The result is shown in Figure [4] as solid lines. The 
best-fit parameters are listed in TableQ] The ionization timescale r 
of on-SNR is small which can explain why CIE model could not 
fit the spectra as well. 


1995) with the extinction Ay = 0.83 magnitude l lStupar et al. 
2007). The elemental abundance is set to the solar value. This 


4 DISCUSSION 

4.1 Distance to SNR G332.5-5.6 

Reliable dis tance to G33 2.5-5.6 is essential to estimate its basic 
parameters. IStupar et alj i2007l) mentioned G332.5-5.6 should 
be within the dust layer of the Galactic disc due to the rich 
shock-excited lines in its spectrum. Therefore it is feasible to build 
an extinction-distance relation to give restrictions to the distance 
of G332.5-5.6. 

Red clump stars (RCs) are core helium-burning early K giants 
which have obvious concentration region in a color-magnitude 


diagram (CMD, e.g. J — K vs. K). Previous studies have 
shown that all RCs nearly have a uniform absolute magnitude 
and their infrared colors are insensitive to metal abundance. In 
view of this, RCs have been used as inf rared stan dard candles 
dLopez-Corredoira et alj200p .lDurant & van Kerkwii l2006l) . 

From the 2MASS All-Sky Point Source Catalog, we extract 
stars surrounding G332.5-5.6 in an area of 0.5 deg 2 with an 
interval of Al = 1.0° and Ab = 0.5°, a total number of 92020 
stars. As seen in Figure[5la), the region mainly occupied by RCs is 
distinct. We divide the K band between 9 and 13 magnitude into 
8 strips with a bandwidth of 0.5 magnitude. For each strip, the star 
counts are calculated in each (J — k) bin of 0.3 magnitude. To 
decrease the influence of unrelated field stars, we only retain 5 bins 
on each side of the bin with maximum star counts. If the bin next 
to the maximum bin has similar counts with the maximum one, we 
will extend 5 bins to 6 bins. Figure [5jb) shows an example with 
the strip K = 11.5-12.0 magnitude. In the example, the equation 
(Durant & van Kerkwiikll2006l) : 


y = A coun t(J - K) a +A RCs exp{——-—~ 

( 1 ) 

is used to fit the star counts distribution. The best-fit parameters 
are (J — K) RCg = 0.81 and a = 0.06. We assume the absolute 
K band magnitude of RCs is Mk = -1.60 to -1.65 mag and 
the in t rinsic color of (J — K) n = 0,63 to 0 . 67 mag (values of 
lAlvesI 200(1 Grocholski & Saraiedinil I2002L lLanev et aD 12012l 


Yaz Gdkce et al. 20j3|). The extinctio n and distance follow 
(Giiyer et al . 2010llL6pez-Corredoirall20T4j) : 


and 


Ay(mag) = 


A k 0.67[(J - K) RCs - (J - K) 0 ] 


0.1137 0.1137 

D{p C ) = lQ^ m K-M K +5-A K )/5] 


( 2 ) 

(3) 


Figure [5jc) shows the extinction-distance relation along 
the line of sight towards G332.5-5.6. The dashed line represents 
the observed averaged extinction of G332.5-5.6. The extinction 
increases rapidly before 3.8 kpc and becomes nearly flat after that 
distance. Since the extinction is mainly caused by the dust, the 
extinction will grow quickly in the dust layer and grow slowly 
outside the dust layer. Therefore, the turnover can be used to 
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Table 1. Best-fit parameters for SNR G332.5-5.6 with model of vnei, vapec and vmekal. a 


vnei 


A r /j(10 22 cm 2 ) 

kT (KeV) 

O (solar) 

Mg (solar) 

Fe (solar) 

T 

Norm (10 4 cm -5 ) 6 

Reduced x 2 /dof c 






(10 10 s cm -3 ) 



0.15 (fixed) 

0.481°;^ 

0 cq+ 0.04 
u - J ^-0.04 

1 70+0.14 
l.^_o.14 

0 76+ 0 - 02 

u./o_o os 

7 C7+O.6I 
z.jj_o 44 

8-00±S;gS 

1.28/533 

rH 

o o 

o o 

+ 1 

o 

0.49±°;°! 

U>:) °-0.05 

1 77+0.14 

0.14 

n 77+O.O6 

U>/ -0.05 

7 C7+O.65 

z - JJ -0.54 

7 . 64 l 3 ;° 3 

1.29/532 

vapec 

Nh(10 22 cm 2 ) 

kT (KeV) 

O (solar) 

Mg (solar) 

Fe (solar) 


Norm (10~ 2 cm -5 ) 

Reduced \ 2 /dof 

0 31+ 0 - 02 
u,J1 -0.03 

O 

OO 

1 + 
o o 

o o 
o o 

0 54+ 0 - 06 

0.05 

7 1 7+0.27 
z - iz — 0.27 

2 30 +0 ' 80 

z.jU—o 58 


1 64 +0 - 25 

0.31 

1.43/533 

vmekal 

NH(l0 22 cm 2 ) 

kT (KeV) 

O (solar) 

Mg (solar) 

Fe (solar) 


Norm (10 -2 cm -5 ) 

Reduced \ 2 /dof 

n 97+O.O2 
U-Z/ -0.02 

n 1 q+ 0-002 
u,l °—0.002 

o 

L/i 

1 + 
o o 

o o 
4^ Cn 

771 +0.28 

1 —0.26 

1 (fixed) 


1 77+O.I7 

A ’ z/ -0.16 

1.53/535 


a Errors are at the 90% confidence level. 

b Norm = ^ D 2 f nenadV, where D is the distance to G332.5-5.6 (cm), n e and n# (cm -3 ) are the electron and hydrogen densities respectively. 
c dof means the degrees of freedom. 



0.0 0.5 1.0 1.5 0.63 0.66 0.69 0.72 0.75 0.78 0.81 0.84 0.87 0.90 0.93 0.96 0.99 1 2 3 4 5 6 7 8 
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Figure 5. (a): CMD for 92020 stars within 0.5 deg 2 of G332.5-5.6. The yellow diamonds and yellow lines show the fitted location of the RCs peak and its 
extent with I a. (b): Histogram of the J — K values of stars extracted from (a) in a 0.5 magnitude interval around K = 11.75 mag. The black curve is the 
best-fit to this histogram using eq. (1). (c): The derived V band extinction-distance relation along the direction towards G332.5-5.6 with overlapped Ay value 
(dashed line) of 332.5-5.6. The red one is for Mjy—- 1.65 mag and (J — k) o=0.67 mag. The black one is for My~- 1.60 mag and (J — fc)o=0.63 mag. 


indicate the boundary of dust layer. Because the observed rich 
shock-excited lines favor the remnant within the dust layer, 3.8 
kpc should be an upper distance limit for the remnant. The lower 
distance limit is directly derived from the relation marked by 
the black points, i.e. 2.2 kpc. Considering all above, we suggest 
a distance of 3.0 ± 0.8 kpc for SNR G332.5-5.6. This distance 
is generally in agreement with previous distances (3-4 kpc, 
iRevnoso & Greenff2007lIstupar et al. 11 20071 ). 


4.2 The Origin of Central X-ray Emission 

Central X-ray emission could be formed by 1) nonthermal X-ray 
emission from a PWN, 2) thermal X-ray emission from the 
shocked ejecta, 3) enhanced central X-ray emission from thermal 
conduction effects which smooth the temperature and electron 
density distribution, 4) projection effect from X-ray emission of 
shocked interstellar medium and 5) evaporation of residual clouds 
inside the remnant. 

IRevnoso & Greenl d2007l) and Istupar et al.l d2007t) had ruled 
out the possibility of nonthermal X-ray emission from a PWN 


because 1) the radio spectral index of the central region is not 
flatter than the index of two outer filaments 2) no pulsar has been 
discovered in the vicinity of G332.5-5.6 3) ROSAT data reveals 
thermal nature of the X-ray emission. By analyzing the Suzaku 
data, we find the central X-ray emission is thermal dominance 
with a subsolar, solar or slightly supersolar elemental abundances. 
This supports neither the PWN origin nor the shocked ejecta origin 
because the heavy elements of ejecta are usually overabundance. 
Thermal conduction, via Coloumb collisions between electrons 
and ions, can reduce the temperature and increase the density 
of the central region of the SNR. This has been suggested to 
explai n the central filled X-ray morphology of mixed-morphology 
SNRs J Shelton et alJl999llSlane et ali200l ). However this model 
usually prod uces a relatively fl at surface brig htness profile (see 
Figure 13 of lShelton et alj[l999l and Figure 6 of Isiane et al.ll2002h 
which can n ot explain the stee p X-ray profile of G332.5-5.6 (see 
Figure 18 of lStupar et al.ll200% . For projection effect and residual 
clouds evaporation, we can not distinguish them because both of 
them are harmonic with the observed data, i.e. element abundances 
and the obvious Ha e mission in the X-ray emission region (Figure 
3 of lStupar et alj 20071) . 









































































Suzaku observation and distance of G332.5-5.6 5 


Table 2. Parameters of SNR G332.5-5.6 


Parameters 

projection Model 

C/t = 2.0 

Evaporation Model 
C/r = 2.5 

C/t = 3.0 

Post-shock Temperature (keV) 

0.49 

0.43 

0.54 

0.68 

Shock radius (pc) 

13.1 

13.1 

13.1 

13.1 

Shock velocity (km s —1 ) 

650 

600 

680 

760 

Age (kyr) 

8 

9 

8 

7 


4.3 The Nature of SNR G332.5-5.6 

Suzaku X-ray data of G332.5-5.6 can be described by plasma with 
temperature of 0.49 keV (5.7 x 10 6 K). This suggests the remnant 
is likely still in the adiabatic phase dKoo et ahll 19951) . 


For the origin of projection effect, a shell SNR evolved in the 
uon-uniform and dense medium can form the mixed morphology 
(this model has tut internal density profiles similar to those in the 
Sedov model, Petruk] |200ll) . Assuming that hydrogen and helium 
have a cosmic abundance ([He]/[H] = 0.1) and all of them are 
ionized, we get n e = 1.2n#. The on-SNR region has an extent of 
about 7\ We assume a depth of l^ (1/12 of the shock radius) and 
a filling factor, /, for the plasma. Under the distance of 3.0 kpc, 
the averaged post-shock electron density will be 0.53// cm ' 3 or 
an averaged total particle density of 0.99/ / cm ~ 3 . There exist the 
Rankine-Hugoniit jump conditions l lRevnoldsll2008l) . 


p2 _ Ui _ 7+1 

pi U2 7—1 


(4) 


rr ~ 2 7(7 - 1) trm p 2 3 pm p 2 

T2 ~ (7+ i)2 Ul= 16 — Ul ' (5) 

here, p, u and T mean density, velocity and temperature respec¬ 
tively. The subscripts 1 and 2 are for per-shock and post-shock. 
For monatomic plasma with solar abundances, 7 = 5/3 and p, = 
0.6, the density of ambient medium is no ~ 0.25// cm ' 3 and 
the shock velocity V s = Ui ~ 650 km s . Taking the radius 
R = 6D = 13.1 pc, we derive a characteristic age t c = I « 8 
kyr. 


For the evaporation model (IWhite & Lond 1 199 ill , if we 
assume the density ratio between the cloud and the intercloud 
medium is large (C —¥ 00 ) and the ratio of the evaporation 
timescale to the age of G33 2.5-5.6 is large ( r —» 00 ), the model 
will only depend on C/t I lKoo et ali 119951) . The central X-ray 
morphology can be formed if C/t > 2. For C/t = 2.0, 2.5 and 
3.0, the post-shock temperature is equal to 0.43, 0.54 and 0.68 
keV which indicate a shock velocity of about 600, 680 and 760 
km s" . The characteristic age of G332.5 will be about 9, 8 and 7 
kyr respectively. We list the derived parameters from both cases in 
Table |T] 


5 SUMMARY 

We have analyzed the Suzaku XIS data of SNR G332.5-5.6. We 
find that the morphology of the central X-ray emission is resembles 
the radio emission. The X-ray data can not be fitted by a CIE 


model, neither vapec nor vmekal, but can be well described by a 
single NEI model, vnei, with a temperature of 0.49jtg Qg keV 
and subsolar abundances for O, Fe and slightly overabundance for 
Mg. The ionization timescale is 2 . 53 jjg '54 x 10 10 cm ' 3 s which 
implies that the plasma is far from ionization equilibrium. This 
explains why CIE model can not fit the X-ray data well. 

Using RCs as standard distance candles, we have derived 
a relationship between extinction and distance along the line of 
sight towards G332.5-5.6 and have estimated its distance of 3.0 ± 
0.8 kpc. We discuss the origin of the central X-ray emission and 
find that the projection effect and evaporation of residual clouds 
likely cause the X-ray emission. Under both projection model and 
evaporation model, we derive the basic parameters of G332.5-5.6 
(Table |2. 
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